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Amendments to the Claims: 

This listing of claims will replace all prior versions, and listings, of claims in the application: 
Listing of Claims: 

1 . (Original) A ligand screening apparatus which screens for a ligand that binds to a protein 
when coordinate data of the protein of single chain or plural chains is given, the apparatus 
comprising: 

post-structural-change protein coordinate data selecting unit that effects structural 
change in consideration of dynamic behavior using induced-fit parameter reflecting induced 
fit on the coordinate data of protein and selects post-structural-change protein coordinate 
data; 

spatial point designating unit that designates a spatial point at which superposition 
with the ligand is to be conducted, from the post-structural-change protein coordinate data 
selected by the post-structural-change protein coordinate data selecting unit; 

interaction function calculating unit that calculates an interaction function when the 
protein and the ligand bind to each other using the spatial point designated by the spatial point 
designating unit and a ligand coordinate data of the ligand; and 

ligand evaluating unit that evaluates the ligand that binds to the protein based on the 
interaction function calculated by the interaction function calculating unit. 

2. (Original) The ligand screening apparatus according to claim 1, wherein the interaction 
function calculating unit calculates the interaction function using Score(iJ) shown in Formula 
1. 



when i is not equal to j 




[Formula 1] 



when i is equal to j 



ax (1-/7) 
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(wherein dy S is a distance between i-th spatial point and j-th spatial point in the target protein. 
dij C is an interatomic distance between i-th atom and j-th atom in the compound, a is a 
coefficient for making Sscore(i j) the maximum value when the group of spatial points in the 
target protein and the compound completely overlap with each other, p is a coefficient for 
giving a threshold value by which it can be defined as "overlapping") 

3. (Currently amended) The ligand screening apparatus according to claim 1 e^2, 
wherein the interaction function calculating unit further comprises interaction function 
optimizing unit that carries out optimization so as to make the score of interaction function 
maximum. 

4. (Original) The ligand screening apparatus according to claim 3, wherein the interaction 
function calculating unit further comprises: 

interaction energy optimizing unit that calculates interaction energy with the protein 
for the superposed ligand after optimization of the interaction function by the interaction 
function optimizing unit, and optimizes the interaction energy while finely adjusting 
conformation of ligand 3D structure data. 

5. (Original) The ligand screening apparatus according to claim 4, wherein the ligand 
evaluating unit further comprises: 

reevaluating unit that executes the interaction function calculating unit after largely 
changing conformation of ligand 3D structure data following optimization by the interaction 
energy optimizing unit, and reevaluates the ligand that binds to the protein based on the 
interaction function calculated by the interaction function calculating unit. 

6. (Currently amended) The ligand screening apparatus according to any on e of claims 
1 te-#, wherein in calculation of any one of the induced-fit parameter and the post-structural- 
change protein coordinate data or both, the post-structural-change protein coordinate data 
selecting unit calculates normal mode for the protein coordinate data, determines intensity of 
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fluctuation of each amino acid, and conduct molecular dynamic calculation using the intensity 
of fluctuation as a constraint condition. 

7. (Original) The ligand screening apparatus according to claim 6, wherein the post- 
structural-change protein coordinate data selecting unit calculates a fluctuation value of 
dihedral angle of main chain according to normal mode calculation, and conducts molecular 
dynamic calculation while setting the fluctuation value as a coefficient of force K in the 
molecular dynamic calculation shown by Formula 2 or Formula 3. 

Erot = Krot (<|> - (|>0) 2 [Formula 2] 

(wherein Erot represents energy of dihedral angle of main chain atom in 3D structure of a 
protein. <() represents dihedral angle of main chain atom. <j)0 represents standard value of 
dihedral angle of main chain atom. Here, when a value of Krot is large, § is constrained by 

*o.) 

Epos = Kpos (r - rO) 2 [ Formula 3 ] 

(wherein Epos represents position energy of main chain atom in 3D structure of a protein, r 
represents coordinate of main chain atom. rO represents standard value of coordinate of main 
chain atom. Here, when a value of Kpos is large, r is constrained by rO.) 

8. (Currently amended) The ligand screening apparatus according to any one of claims 
1 te-^, wherein the interaction function calculating unit uses the interaction function to which 
a dynamic property function representing dynamic property of protein is added as "elastic 
energy". 

9. (Original) The ligand screening apparatus according to claim 8, wherein the interaction 
function calculating unit adapts "U collision" as "elastic energy" which is a function shown 
by Formula 4 in consideration of local flexibility of protein. 
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M N 

U collisi on =ZS^( / '>) 

<p (/, y ) = Kcollision * (Rcollision (i, j) -R) 2 

[Formula 4] 

(wherein M is number of atoms in an active site that prohibit collision, N is number of atoms 
of ligand. When interatomic distance R between an i-th atom of a main chain or a side chain 
with a little dynamic behavior in an active site, and j-th atom in the ligand is not more than 
collision distance "Rcollision(ij)", (|>(i j) is calculated) 

1 0. (Currently amended) The ligand screening apparatus according to any on e of claims 
1 te-7, wherein the interaction function calculating unit uses the interaction function to which 
a normal mode analysis result or secondary structure determination result of the protein is 
added as a dynamic property function that represents dynamic property of protein. 

1 1 . (Original) A ligand screening method which screens for a ligand that binds to a protein 
when coordinate data of the protein of single chain or plural chains is given, the method 
comprising: 

post-structural-change protein coordinate data selecting step that effects structural 
change in consideration of dynamic behavior using induced-fit parameter reflecting induced 
fit on the coordinate data of protein and selects post-structural-change protein coordinate 
data; 

spatial point designating step that designates a spatial point at which superposition 
with the ligand is to be conducted, from the post-structural-change protein coordinate data 
selected by the post-structural-change protein coordinate data selecting step; 

interaction function calculating step that calculates an interaction function when the 
protein and the ligand bind to each other using the spatial point designated by the spatial point 
designating step and a ligand coordinate data of the ligand; and 



WASH_1671616.1 



-5- 



Atty. Dkt. No. 079088-0102 



ligand evaluating step that evaluates the ligand that binds to the protein based on the 
interaction function calculated by the interaction function calculating step. 

12. (Original) The ligand screening method according to claim 11, wherein the interaction 
function calculating step calculates the interaction function using Score(i j) shown in Formula 
1. 



(wherein d x f is a distance between i-th spatial point and j-th spatial point in the target protein. 
djj C is an interatomic distance between i-th atom and j-th atom in the compound, a is a 
coefficient for making Sscore(i j) the maximum value when the group of spatial points in the 
target protein and the compound completely overlap with each other, p is a coefficient for 
giving a threshold value by which it can be defined as "overlapping") 

13. (Currently amended) The ligand screening method according to claim 1 1 or 12 , 
wherein the interaction function calculating step further comprises interaction function 
optimizing step that carries out optimization so as to make the score of interaction function 
maximum. 

14. (Original) The ligand screening method according to claim 13, wherein the interaction 
function calculating step further comprises: 

interaction energy optimizing step that calculates interaction energy with the protein 
for the superposed ligand after optimization of the interaction function by the interaction 
function optimizing step, and optimizes the interaction energy while finely adjusting 
conformation of ligand 3D structure data. 



when i is not equal to j 




[Formula 1] 



when i is equal to j 
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15. (Original) The ligand screening method according to claim 14, wherein the ligand 
evaluating step further comprises: 

reevaluating step that executes the interaction function calculating step after largely 
changing conformation of ligand 3D structure data following optimization by the interaction 
energy optimizing step, and reevaluates the ligand that binds to the protein based on the 
interaction function calculated by the interaction function calculating step. 

16. (Currently amended) The ligand screening method according to any one of claims 1 1 
to 15 , wherein in calculation of any one of the induced- fit parameter and the post-structural- 
change protein coordinate data or both, the post-structural-change protein coordinate data 
selecting step calculates normal mode for the protein coordinate data, determines intensity of 
fluctuation of each amino acid, and conduct molecular dynamic calculation using the intensity 
of fluctuation as a constraint condition. 

17. (Original) The ligand screening method according to claim 16, wherein the post- 
structural-change protein coordinate data selecting step calculates a fluctuation value of 
dihedral angle of main chain according to normal mode calculation, and conducts molecular 
dynamic calculation while setting the fluctuation value as a coefficient of force K in the 
molecular dynamic calculation shown by Formula 2 or Formula 3. 



(wherein Erot represents energy of dihedral angle of main chain atom in 3D structure of a 
protein. § represents dihedral angle of main chain atom. <|>0 represents standard value of 
dihedral angle of main chain atom. Here, when a value of Krot is large, $ is constrained by 




[Formula 2] 



<j>0.) 




[Formula 3] 
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(wherein Epos represents position energy of main chain atom in 3D structure of a protein, r 
represents coordinate of main chain atom. rO represents standard value of coordinate of main 
chain atom. Here, when a value of Kpos is large, r is constrained by rO.) 

1 8. (Currently amended) The ligand screening method according to any on e of claims 1 1 
to 17 , wherein the interaction function calculating step uses the interaction function to which 
a dynamic property function representing dynamic property of protein is added as "elastic 
energy". 

19. (Original) The ligand screening method according to claim 18, wherein the interaction 
function calculating step adapts "U collision" as "elastic energy" which is a function shown 
by Formula 4 in consideration of local flexibility of protein. 

M N 

U coms ion =E2>0W) 

<p (/, j) = Kcollision * (Rcollision (i, j) -R) 2 

[Formula 4] 

(wherein M is number of atoms in an active site that prohibit collision, N is number of atoms 
of ligand. When interatomic distance R between an i-th atom of a main chain or a side chain 
with a little dynamic behavior in an active site, and j-th atom in the ligand is not more than 
collision distance "Rcollision(ij)", (|>(i j) is calculated) 

20. (Currently amended) The ligand screening method according to any one of claims 1 1 
to 17 , wherein the interaction function calculating step uses the interaction function to which 
a normal mode analysis result or secondary structure determination result of the protein is 
added as a dynamic property function that represents dynamic property of protein. 

21 . (Original) A program which makes a computer execute a ligand screening method which 
screens for a ligand that binds to a protein when coordinate data of the protein of single chain 
or plural chains is given, the method comprising: 
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post-structural-change protein coordinate data selecting step that effects structural 
change in consideration of dynamic behavior using induced-fit parameter reflecting induced 
fit on the coordinate data of protein and selects post-structural-change protein coordinate 
data; 

spatial point designating step that designates a spatial point at which superposition 
with the ligand is to be conducted, from the post-structural-change protein coordinate data 
selected by the post-structural-change protein coordinate data selecting step; 

interaction function calculating step that calculates an interaction function when the 
protein and the ligand bind to each other using the spatial point designated by the spatial point 
designating step and a ligand coordinate data of the ligand; and 

ligand evaluating step that evaluates the ligand that binds to the protein based on the 
interaction function calculated by the interaction function calculating step. 



22. (Original) The program according to claim 21, wherein the interaction function 
calculating step calculates the interaction function using Score(ij) shown in Formula 1. 



Sscore(i,j) = ^ 



when i is not equal to j 
ax[exp {-{d;-dff}-j3_ 



when i is equal to j 
ax(l-j3) 



[Formula 1] 



(wherein dy S is a distance between i-th spatial point and j-th spatial point in the target protein, 
dy is an interatomic distance between i-th atom and j-th atom in the compound, a is a 
coefficient for making Sscore(i j) the maximum value when the group of spatial points in the 
target protein and the compound completely overlap with each other. P is a coefficient for 
giving a threshold value by which it can be defined as "overlapping") 



WASH_1 671 616.1 



-9- 



Atty. Dkt. No. 079088-0102 



23. (Currently amended) The program according to claim 21 or 22, wherein the 
interaction function calculating step further comprises interaction function optimizing step 
that carries out optimization so as to make the score of interaction function maximum. 

24. (Original) The program according to claim 23, wherein the interaction function 
calculating step further comprises: 

interaction energy optimizing step that calculates interaction energy with the protein 
for the superposed ligand after optimization of the interaction function by the interaction 
function optimizing step, and optimizes the interaction energy while finely adjusting 
conformation of ligand 3D structure data. 

25. (Original) The program according to claim 24, wherein the ligand evaluating step further 
comprises: 

reevaluating step that executes the interaction function calculating step after largely 
changing conformation of ligand 3D structure data following optimization by the interaction 
energy optimizing step, and reevaluates the ligand that binds to the protein based on the 
interaction function calculated by the interaction function calculating step. 

26. (Currently amended) The program according to any on e of claims 21 to 25 , wherein 
in calculation of any one of the induced-fit parameter and the post-structural-change protein 
coordinate data or both, the post-structural-change protein coordinate data selecting step 
calculates normal mode for the protein coordinate data, determines intensity of fluctuation of 
each amino acid, and conduct molecular dynamic calculation using the intensity of fluctuation 
as a constraint condition. 

27. (Original) The program according to claim 26, wherein the post-structural-change protein 
coordinate data selecting step calculates a fluctuation value of dihedral angle of main chain 
according to normal mode calculation, and conducts molecular dynamic calculation while 
setting the fluctuation value as a coefficient of force K in the molecular dynamic calculation 
shown by Formula 2 or Formula 3. 
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Erot = Krot(<|>-<J>0) 2 



[Formula 2] 



(wherein Erot represents energy of dihedral angle of main chain atom in 3D structure of a 
protein. <|) represents dihedral angle of main chain atom. <|>0 represents standard value of 
dihedral angle of main chain atom. Here, when a value of Krot is large, (J) is constrained by 



(wherein Epos represents position energy of main chain atom in 3D structure of a protein, r 
represents coordinate of main chain atom. rO represents standard value of coordinate of main 
chain atom. Here, when a value of Kpos is large, r is constrained by rO.) 

28. (Currently amended) The program according to any on e of claims 21 to 27 , wherein 
the interaction function calculating step uses the interaction function to which a dynamic 
property function representing dynamic property of protein is added as "elastic energy". 

29. (Original) The program according to claim 28, wherein the interaction function 
calculating step adapts "U collision" as "elastic energy" which is a function shown by 
Formula 4 in consideration of local flexibility of protein. 



(wherein M is number of atoms in an active site that prohibit collision, N is number of atoms 
of ligand. When interatomic distance R between an i-th atom of a main chain or a side chain 
with a little dynamic behavior in active site, and j-th atom in the ligand is not more than 
collision distance "Rcollision(ij)", (|)(ij) is calculated) 



+0-) 




[Formula 3] 



U coiiision =Z2M / >- / *) 



/=1 y=l 

<p (/, j) = Kcollision * (Rcollision (/ , j) - R) 2 



[Formula 4] 
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30. (Currently amended) The program according to any one of claims 21 to 27 , wherein 
the interaction function calculating step uses the interaction function to which a normal mode 
analysis result or secondary structure determination result of the protein is added as a dynamic 
property function that represents dynamic property of protein. 

3 1 . (Currently amended) A computer readable recording medium in which the program 
according to e ither of claims 21 to 30 is recorded. 
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